Ro-vibrational CO Detected in the j3 Pictoris Circumstellar Disk 
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ABSTRACT 

We present high resolution near-infrared spectra of f3 Pictoris - a nearby young star with 
a debris disk. Fundamental low-J CO absorption lines are detected and strict upper limits are 
placed on the flux of v=2-l low-J CO emission lines. The limit on the UV fluorescent emission flux 
in the v=2-l lines is used to place a tight constraint on the inner extent of the CO gas. Assuming 
H I is the primary collision partner, the subthcrmal population of the low-J v=0 rotational levels 
constrains the density of the gas in the disk to nn = (2.5ti'3>) x 10 5 cm~ 3 . If the distribution 
of hydrogen follows that of the other metals in the disk, we find that the mass of the gas in the 
disk is (Q.ntoil) M e- We compare this mass to the gas mass necessary to brake the metals in 
the disk through ion-neutral reactions. 

Subject headings: circumstellar matter-infrared: stars, planetary systems: protoplanetary disks 



1. INTRODUCTION 



B Pictoris (f3 Pic) is a 12 Myr (jZuckerman et al.ll200ll ). A5V star at a distance of 19.3 pc (jCrifo et al 



19971 ) with a debris disk. The d ust component of the disk has been studied e xtensively from the op- 
tical through the sub- millimeter (Smith fc Terrile 1984 : Artvmowicz et al. 19891: Golimowski et al. 1993 : 



Kalas fe Jewittlll995tlMouillet et al.ll997t|Pantin etal.ll997tlHolland.et al.ll998tlHeap et alj200 c ULiseau et al 



2003 



Okamoto et al.ll2004i I Golimowski et al 



2006 



bution, and in particular warps in the disk (e.g. lKalas fc Jewitt 



Boccaletti et al 



;t al l 
199.4 



2009) . Asymmetries in the dust distr i 



Mouillet et al 



point to the presence of a planet in the disk which was recently imaged (jLagrange et al. 2009, 201 



1997tlHeap et al.l l2000h. 
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Ferlet et al. 


1987; 


Jollv et al. 


1998; 


Roberee et al. 


2000; 


— . 

Brandeker et al. 


2004) 



red-shifted rela tive to the rest frame of Pic, resulting from infalling gas released from star-grazing plan- 



etesimals (e.g., IVidal-Madiar et al.l ll994). A second component (comprised of Na I, Fe I, Ca II, Ti II, Ni I, 



Ni II, Cr II) rema i ns at rest relative to the star and has a sp atial distribution similar to that of the dust 



(jFerlet et al.lll987t lOlofsson et alj|2001t iBrandeker et al.l 12004'). The existence of the stable gas in a debris 



disk was somewhat surprising because radiation p ressure from Pic to should be sufficient to radially ac- 



celerate the gas to high velocities (> 100 km s 1 ; IBrandeker et al.ll2004f). Thus, some braking mechanism 



is necessary to explain the stable component of the gas. IBrandeker et al. ( 2004 ) e stima ted that ~ 50M® 
of gaseous hydrogen is needed if it is the braking agent. However. iFreudling et al.l (|1995f ) show from limits 
on 21 cm emission that the atom ic hydrogen mass must be less than 1.6 Mm. As for the molecular gas, 
Lecavelier des Etangs et al.1 (|200lh place a tight constraint on the mas s of H2 at < -l Mm based on the non- 
detection of absorption lines from electronic transitions. Additionally, IChen et al.l (|2007T ) use the absence of 
H2 emission in Spitzer IRS spectra to limit the a mount of warm H2 (50 - 1 00 K) in the disk at <17 M®, ruling 
out a more extended reservoir of molecular gas. iThebault fe Augereaul (|2005h use dynamical arguments to 
argue that the total gas mass must be <0.4 M®. Although the se observational c onstra ints would appear 
to suggest that the disk is not massive enough to brake the gas, IFernandez et al.l ((2006J) note that a much 
smaller mass of hydrogen may be sufficient to brake the gas because only the ionized component needs to be 
decelerated. This is because the elements that receive the highest radiative force are also the elements with 
the highest ionization rates, th us only the ions are imp ortant. The acceleration of the neutrals is negligible, 
thus they do not need braking. IFernandez et al.l ( 20061 ) explored possibilities for braking the ions in the gas, 



which would explain the observed stable component of the gas. 

The first possibility is ionized particle collisions, where the ions lose momentum by Coulomb interactions. 
The ions are dynamically coupled, so the radiative force on the particles can be described by an average over 
all of the particles. Some ions feel a low radiative force, so if these atoms were somehow overabundant, the 
average radiative force on the ions would go down - effectively slowing the gas. Carbon is one such atom 
that feels a low radiative force. Indeed, an increase in the carbon abundance by a factor of ~10 with respect 
to other atoms may be sufficient to lower the average radiative force on the ions, braking the gas. Atomic 
carbon is observed to be enhanced compared to oxygen in the Pic system by a factor of 18 relative to solar 
|Roberge et al.ll2006l h making this an attractive possibility. 



Another possibility for braking the gas is ion collisions w ith neutral gas. In this scenario, a total 
disk mass of at least 0.1 Mm is required ( Fernandez et al. 20061) . well within all of t he ob servational limits 
mentioned above. This is much smaller than the value inferred by IBrandeker et al.l (|2004f ). as they did not 
assume only the ionized particles need to be slowed. While the ion-ion collision scenario is possible, limits on 
the disk mass cannot rule out a contribution to the stabilization of the atomic gas from ion-neutral collisions 
as well. 

To explore the role ion-neutral collisions play in braking the gas in the outer disk around Pic, we present 
near infrared high resolution spectra of Pic centered near 4.7 /xm. We detect fundamental ro-vibrational 
CO lines in absorption and place strict constraints on fluorescent emission from higher vibrational bands. We 
find that the rotational levels are only thermalized up to J"=2 and describe how this information provides a 
means to calculate the density of gas in the disk. We also use the non-detection of fluoresced ro-vibrational 
CO emission to constrain the inward extent of CO in the disk. 
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OBSERVATIONS 



W e acquired high-resolution near-infrared spectra of Pic using PHOE NIX dHinkle et al.ll2003l l2000l 
1998) at the Gemini South telescope and CSHELL (jTokunaga et al.l ll99oh at the NASA Infrared Tele- 



scope Facility. The resolutions of PHOENIX and CSHELL are R=50,000 and R~40,000, respectively. The 
PHOENIX observations were taken on March 23, 2008 and January 12, 2003. The CSHELL observations 
were taken on August 8, 2000. The observations were centered around 4.7 fim to cover fundamental ro- 
vibrational transitions of CO. A summary of observations, presented in Table [TJ includes the lines that are 
observed. 

Observations in the 4.7 /xm region are dominated by a strong thermal background, limiting the exposure 
times. Therefore, short exposures are taken while nodding between two positions separated by ~5". The 
exposures are taken in an ABBA pattern in order to cancel the thermal continuum to first order. The 
scans are flat fielded, cleaned of hot and dead pixels, including cosmic ray hits, and then combined in the 
sequence (Ai-Bi-B2+A2)/2. Because the spectra are curved along the detector, they are first rectified by 
finding the centroid of each column and shifted to a common row. A one-dimensional spectrum is then 
extracted from the rectified spectrum. This spectrum is combined with an atmospheric transmittance model 
spectrum in order to find a wavelength solution. The model is created using the Spectral Synthesis P rogram 
|Kunde k Maguirelll974l ). which accesses the 2000HITRAN molecular database (jRothman et alj|2003l) . Each 
spectrum is then ratioed to a standard star observed at a similar airmass to remove telluric absorption lines. 
Areas where the transmittance is below 50% are omitted. To determine the signal to noise of our spectra, 
we measured the standard deviation of continuum (Table Q]). The extracted spectra are presented in Figures 
[T] and [5J The signal to noise of the spectrum acquired with PHOENIX was 100. The spectrum acquired 
with CSHELL is dominated by the broad Pf feature and mis-canceled water line. The signal to noise of 
the spectrum acquired with CSHELL was 90. 



RESULTS 



We detect the R(0), R(l), and R(2) fundamental ro-vibrational CO absorption lines near 4.64 /jm 
(Figured]). We do not detect the R(3) feature. The R(l) and R(2) lines are unresolved at the the 6 km 
s _1 resolution of PHOENIX, which is consistent with the intrinsic line width b=1.3 km s" 1 inferred from 



the electronic absorption measurements observed in the UV (jRoberge et al.l 120001 ) . Similarly, the R(0) line 
observed with CSHELL was unresolved at the effective ~8 km s _1 resolution of the instrument. The R(l) 
was not detected with CSHELL as the CO line was not cleanly separated from the heavily saturated telluric 
component. The N2O line on the red wing of the telluric feature lowers the transmittance of the atmosphere 
an additional 20%. The heliocentric Doppler shift of the CO absorption lines for all dates arc +21 ± 1 
km s" 1 , consistent with the heliocentric Doppler shift of Pic, and indicates that it is part of the stable 



component of the gas disk (Ferlet et al. 1987 ; Jolly et al. 1998 ; Roberge et al. 2000l : Brandeker et al. 2004 : 



Roberge et al.l 2006; Roberge k Wci nberge: 



200 



Figure [2] shows the spectrum in the 4.74 /im region containing the v=l-0 P(8) and P(9) lines and the 
v=2-l P(l), P(2), and P(3) lines. These features are not detected. While the gas in the disk around Pic 
is too cold to give rise to collisionally excited CO ro-vibrational e mission or absorption, CO ro-vibrational 
emission can be excited by UV fluorescence ([Krotkov et al.lll980l ). When CO is excited electronically and 
relaxes back to the ground electronic state, the excited vibrational bands are populated. Thus, even relatively 
cool gas can give rise to a large vibrational temperature. For example, the CO around HD 141569 has a 
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rotational temperature of 200 K while the vibrational temperature is ^5000 K ( Brittain et al. 20031 2007 ). 
CO was first detected around j3 Pic in the UV via electronic absorption lines. Thus, one may expect to 
observe cold ro- vibrational emission from /? Pic. We do not detect this emission in our spectrum around 4.74 
fim, with a signal-to-noise of 90. There appears to be a spurious peak at the at the location of the v=l-0 
P(9) feature. However, this is not assumed to be a true featue. The peak is embedded in a combination of 
telluric features, with the transmittance at the location only at 68%. Additionally, the P(9) emission cannot 
be significantly brighter than the non-detected P(8) feature. Finally, the gas would have to be quite warm 
to detect it this far up the rotational ladder, which is also not observed (see Section [472]) . Assuming the lines 
are Gaussian and unresolved, our one-sigma upper limit on the equivalent width (based on the S/N of the 
spectrum) is ~ 5 x 10~ 4 cm -1 . Scaling the equivalent width by the flux density as measured by the Infrared 
Space Observatory (F4.7 Mm = 1 x 10~ 9 ergs s _1 cm~ 2 /im -1 ), the upper limit on the flux of the individual 
lines is 7.1 x 10" 16 ergs cm" 2 s" 1 (Table©. 



ANALYSIS 



4.1. Constraint on CO Distribution 



The excitation of the electronic transitions of CO via the absorp tion of UV flux resu lts in the popu- 
lation of excited vib rational levels in the ground electronic state (e.g. iKrotkov et al.lll980l ). For example, 
Brittain et al.l (|2003l) observe ro-vibrational CO emission from the disk around HD 141569 distributed from 
the inner 50 AU and conclude that the vibrational temperature of the gas is ^5500 K while the rotational 
temperature of the gas is ^200 K. This situation occurs when the vibrational population of cool gas is domi- 
nated by the de-excitation of electronically excited molecules. In this "strong pumping limit," the vibrational 



popu lation of CO near the star reflects the color temperature of the star's ultraviolet field ([Krotkov et al 



1980). Further from the star, the spontaneous de-excitation of the vibrational levels becomes more significant 
relative to the rate at which these levels are populated by the relaxation of electronically excited molecules 
and the vibrational temperature decreases. 

We can use the non-detection of the CO emission in the disk to constrain the distribution of CO using 
the upper limits on the v=2-l emission features (displayed in Figure [5]) and the analysis in the previous 
section. The fluxes of the v=2-l P(l), P(2) and P(3) lines are < 7.1 x 10~ 16 erg cm" 2 s" 1 (Section[3|). 



Jolly et al.1 (|1998l ) and lRoberge et al.l (|2000l ) observe the CO absorption lines near 1500A that represent 
transitions from the ground electronic state to the first excited electronic state. When these molecules relax, 
they will populate v>l in the ground electronic state. We have calculat ed the population of v =2 for CO as 
a function of distance from the star following the procedure described by Brittain et al. (2009). We take our 
input UV spectrum for (3 Pi c taken from Kuru cz atmospheric models (Kurucz 19931 ) . which are then scaled 
to the emission presented in I Jolly et al. (1998). 



As a limiting case, we fi rst assume that the CO in the disk is distributed over the same radial range as 
the atomic gas observed by (jBrandeker et al.l 12004). The UV is absor bed over an annulus extending from 
13 AU to 323 AU, with a profile determined bv lBrandeker et al.1 (|2004l ). They infer the radial profile for the 
gas disk based on their measurement of the distribution of Na I (discussed in detail in Section 14. 3[) : 

-1/2 



nco(r) = n 



/ r / r \ 

V 117 AU/ + \ 117 AU/ 



5.3 



(1) 



where n is the fiducial density of CO at 117 AU. 
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The UV flux propagates through each gas annulus, where the flux is attenuated by the distance from 
the host star as it travels through the annulus. Throughout the annulus, the population of CO molecules in 
v=2 falls of with the distance to the star. In this case, no is normalized to reproduce a line-of-sight column 
density of 2.1 x 10 15 cm -2 (calculation of the column density is described in Section B~2"1) . We require an 
inner radius rj n > 25 AU for the disk to intercept a small enough fraction of the stellar flux that the emergent 
UV fluorescent flux is below the observational limits. If the CO is restricted to a smaller radial extent, the 
number of molecules in v=2 will increase as more gas is nearer the star. In this case, the inner edge of the 
gas is at a larger radius. Thus, the inner radius must be larger than 25 AU. The assumption that the CO 
has the same radial distribution as the Na I is a limiting case. If a different profile is used, for example 
something closer to a Kuiper Belt distribution where the density of CO would fall off faster than the profile 
used, then the CO population in v=2 would increase closer to the star, which would increase the limit on 
the inner radius of CO emission to remain consistent with the non-detection of the v=2 features. Thus, the 
inner radius has a hard constraint at ri n > 25 AU. There is no constraint on the outer radius (Figure [3]). 



4.2. Ro-vibrational Absorption 

We can also use the equivalent width of the absorption features to determine the rotational temperature 
of the absorbing gas. The column density of molecules in each state is related to the equivalent width of 
absorption by: 

Nj " = 8.85x10-"/^,, (2) 
where J" is the lower transition state, Wp is the equivalent width of the line, 8.85 x 10 -13 cm is the classical 
electron radius (7re 2 /m e c 2 ), and fj>j>> is the oscillator strength of the transition. The equivalent widths and 
column densities are presented in Table [2J If the CO is in local thermodynamic equilibrium, then the relative 
population of each level is given by the Boltzmann distribution, 

J\f jr , _ Atot(2./" + 1) c -hcBj"(j" + l)/kT ^ 

where Q is the rotational partition function, and B is the rotational constant. This can be rewritten as: 



hcB V 2 J" + 1 / T ' hcB V Q 

where A to t is the total number of CO absorbing molecules. Thus, plotting (k/hcB) Zn(Aj,,/2J"+l) vs 
J"(J"+1) results in a linear relationship where the negative reciprocal of the slope is equal to the rotational 
temperature. However, when we plot our data (Figure U]), we find that while the first three levels are 
consistent with a rotational temperature of 15 ± 2 K, the column density of molecules in J" = 3 is under- 
populated by a factor of 2.6. A sub-thermal population can occur when the gas density is too low for 
collisions to balance spontaneous emission. The near-infrared transitions are optically thin, so we sum the 
column densities from J" = to J" = 2 to find the column density of CO, (2.1 ± 0.3) x 10 15 cm -2 , with 
the column density of CO in J" = 3 is < 2.3 x 10 14 cm -2 . We assume the overall CO column density to be 
(2.1 ±0.3) x 10 15 cm" 2 . 



Jolly et al.l (|1998l ) used electronic absorption spectra of CO from the Goddard High-Resolution Spectro- 
graph (^1500A) to meas ure the column densit y and temperature of CO and found N(CO) = (2 ± 1) x 10 15 
cm- 2 and T=20 ± 5 K. Roberge et al. ( 2000[) used the Space Telescope Imaging Spectrograph (STIS) to 



find a CO column density of N(CO) = (6.3 ± 0.3) x 10 14 cm" 2 at a temperature of 15.8 ± 0.6 K. Our 
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column density is con sistent with that of iJollv et al.l (|1998l) . and a factor of ~3 larger than that reported by 
Roberge et al. ( 2000l). whi l e the temperature we derive is consistent with both studies. As both Jolly et al 



(|1998l ) and lRoberge et al.l ([20001 ) are measuring absorption of the same electronic transitions, these two UV 
studies should be tracing the same column of gas. Further, the ro-vibrational absorption we observe should 
probe the c olumn of gas obser ved against the continuum of the star as there is no extended Lyct emission 
in the UV (jMarilli et al.lll9971 ). so it is u nclear why the col umn density observed on t hese three epochs is 
different, though the discrepancy betw een Jolly et al. J 19981 ) and Roberge et al. ( 2000 ) is only 1.5a due to 



the large error on the measurement bv lJollv et al 



(199 



One possibili ty is that the abs orption is variable on t he 3 year timescale b etween the UV observations on 
November 1994 ([Jolly et al.lll998l ). and December 1997 (jRoberge et al.ll2000h . We do not see this degree of 
variability between the August 2000 CSHELL and the January 2003 Phoenix observations. If the column of 
CO changed by a factor of three during this time, we would not find a linear rotational temperature diagram. 
Indeed, a change in the column density by a factor of three between the 2000 and 2003 observations would 
cause the J" = point to be offset by ~ 0.4 units relative to the other data points in our excitation 
diagram (FigureHJ). Another explanation may be that th e UV observations of the J" < 3 lines are saturated, 
providing an underestimate to the total column density. Roberge et al. ( 2000l) assumed the rotational levels 
were in LTE; however, our results indicate that this is not the case (see Figure 2]). If the lowest three 
lines were slightly saturated (to > 1), then it is possible that the total column of CO could be significantly 
underestimated by assuming a thermal population while fitting to optically thin, sub-thermally populated 
lines. To determine if the assumption of NLTE can reproduce the observed equivalent widths of the CO 
fundamental lines, we calculate the rotational levels explicitly in the next section. 



4.3. Excitation 

The population of the rotational levels of CO are calculated explicitly under the assumption that the 
molecules are excited by a combination of collisional excitation and radiative pumping by the cosmic mi- 
crowave bac kground, and that th e level populations are in steady state. The lifetime of CO in the disk is 
^200 years (jRoberge et al.1 120001 ) based on the timescale of dissociation from interstellar UV photons, so 
we assume that the creation and destruction of CO in the disk is unimportant. To highlight the role these 
various processes play in determining the populations of each rotational level, we consider the state J=l, 



dn 



J=i 



dt 



-BqiJoiuq — (Aio + B w J w )ni 

7 

- ^ k u n x + (A 2 i + B 2 iJ2i)n 2 

6 

- k ml n m = 0, 



(5) 



where By is the Einstein coefficient, Jij is the energy density of the cosmic microwave background at the 
frequency of the transition, and k^ is the collisional (de)excitation coefficient from the ith state to the jth 
state, and rii is the fractional population of the i th level. The collisional coefficients are given by 



(6) 
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for an upward transition (i < j), and 



t - 2J ' 

01 ~ 2J 



+l e - h cZ /k TJ2n s f 8 (T) 



(7) 



for a downward transition (i > j). 

The variable n s is the number density of the collision partner, J" is the lower state, J' is the upper state, 
and f s (T) is the collision rate for a specific species, s. A suite of synthetic spectra were created with the 
num ber density of the col lision partner as a free parameter. We assume an intrinsic line broadening of 1.3 km 
s- 1 ([Roberge et al.ll2000l) . The equivalent widths of the features from the synthetic spectrum are compared 
to the observed equivalent widths. We determine the best fit to the density by calculating chi-squared for 
each density of a given collision partner. 

The CO in the disk is likely being photodesorbed from icy grains. This is a process by which the 
molecules in the ice lattice are electronically excited. This energy breaks the bonds of the lattice allowing 
the mol ecules to escape. T he photodesorption efficiency of CO is similar to the photodesorption efficiency 
of HoO dOberg et al.ll2009allbh . The mixing ratio for H 2 0/CO can range from as l ow as theoretical values o f 



1 (D'H endecourt et al 



1985 ) to as high as ^50 as observed in Oort cloud comets ( Disanti fc Mumma 20081) 



Thus it is plausible that H2O could be the dominant collision partner with CO. 



The collision rate between H2O and CO has been estimated bv lGreenl (|l993l ) to be 1.04 x 10~ 
at 15 K for the transition between the ground and first excited rotational level. Adopting these parameters, 
we find that the best fitting H2O density is nn 2 o = (3.0^'q) x 10 3 cm 3 - Figure [S] shows the chi-squared 
values assuming 15 K for different densities of H2O. If we assume that the column density of H2O is 1-50 
times the column density of CO, and the CO and H2O are co-spatial in an annulus, then the path length of 
the water annulus would be: 

_ N{H 2 Q) 

-^annulus — \y ) 

nn 2 o 

This leads to a radial path length of the gas annulus ranging from 4.7 x 10 -2 AU to a maximum of 2.3 
AU. The desorption efficiency depends on the intensity of the incident UV, which decreases as r~ 2 in the 
radial direction. It seems unlikely that the desorption would decrease sharply at a distance of 25 AU from 
the star, restricting the gas to such a small annulus. Both CO and H2O should be coming off of grains just 
as prodigiously at, for example, 35 AU, as at the inner radius of 25 AU. While it is true that if the disk 
is extended, it is implausible that H2O is the dominant coll ision partner, observations show evid ence for 



disks that have sculpted, narrow annuli (such as Fomalhaut: iKalas et alJl2008t Ichiang et alJl2009h . Othe 



species, such as CH4, may also photode sorb from grains an d collisionally excite CO, as the photodesorption 
efficiencies are similar to CO and H 2 ( Oberg et al. 2009al lbh . Collision rates between these species and CO 
have not been measured, though if we assume that collision rates for species like CH4 behave similarly to 
molecules such as H 2 and N 2 (the rates for H 2 were calculated from N 2 collisions), the same conclusion 
can be drawn about other molecules being a dominant collision partner. Thus, we cannot rule out the 
possibility that H2O (or other similar molecules) can explain the excitation of CO. 

An alternative scenario is that at an age of of 12 Myr, /3 Pic retains some primordial hydrogen. Thus the 
CO could be photodesorbed from icy grains and collisionally excited by primordial hydrogen gas in the disk. 
While H I has not been detec ted, the upper limits do not rule o ut the possibility that th e hydrogen abundance 
reflects a solar composition (jBrandeker et al] 120041 ) . Indeed. iBrandeker et al. (120041) use a photoionization 
code to infer the radial profile of hydrogen in the gas disk based on their measurement of the distribution of 
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Nal: 



n H (r,z) = no 



/ r \ 2A / r \ 
U17AU/ + V117 AU/ 



5.3 



-1/2 



x exp 



(0.17r) 2 J ' 



(9) 



where the metals are observed to extend from 13 AU to 323 AU. Assuming the disk reflects a solar compo- 
sition, they find the mas s of the disk is ~0-l M m. This is below the detection limits determined from 21 cm 
observations of the disk (jFreudling et al.lll995l ). 



If hydrogen is the dominant collision partner, we can determine the density of H I necessary to col- 
l isionally excite the o bserved rotational levels of CO. We adopt the distrib ution of H I sugges t ed by 
( Brandeker et al. 2004 ) and the collision rates between hydrogen and CO from Green fc Thaddeusl ( 1976j) 
using Equation [5] We assume that the CO originates in an ann ulus embedded in a m uch larger radial dis- 
tribution that includes hydrogen and other metals observed by (jBrandeker et al.ll2004l) . That is, the CO is 
coming off of grains into a bath of diffuse gas, including hydrogen. We do not assume the CO and hydrogen 
have the same inner radius, but instead determine the density of hydrogen at the inner radius of CO (25 AU; 
from Section l4~Tj) . Of course, the CO may originate at a slightly larger radius, and it isn't clear how wide 
the annulus is. The best-fitting density for the hydrogen at 25 AU in order to reproduce the CO populations 
observed is nn = (2. 5^2) x 10 5 cm" 3 (Figure [BJ. Interestingly, the upper limit of the hydrogen density 
at the 3cr level is unconstrained. Thus, assuming the levels are not in LTE sufficiently can describe the 
observations, possibl y resolving the factor of ^3 difference between the CO column density determined by 
Roberge et al] (|2000i ) and the column density presented in this work (Section FT~2"]) . The fit to the absorption 
features is shown in Figure[7J There is a prediction for the R(3) feature, which seems to fit the 1.5er feature 
at the same location. If the feature is determined to be real with further observations, it would more strongly 
constrain our results, but it would not appreciably alter our conclusions. That is, the population would still 
prove to be sub-thermalized. 



4.4. Constraint on Disk Column Density and Mass 



Continuing with this scenario, we assume the CO is embedded in a much larger gas disk (Figure [3]) made 
up of hydrogen and other metals, and that the density of hydrogen we determined is the fiducial density 
at the mid-plane (z=0), and the hydrogen continues inward beyond the CO (>25 AU) to where the atomic 
ga s originates (13 AU). If the hydrogen follows the density profile described using the distribution given 
by iBrandeker et al.1 (|2004l Equation [9]) for other metals in the disk, then the fiducial H I density scaled 
to 117 AU becomes no=(3.9l^ 1 8 ) x 10 4 cm -3 . Note that the hydrogen density calculated in Section l4~3l is 
independent of the radial distribution of CO, but rather is calculated by fitting to the level populations. 
We can calculate the line-of-sight column density of atomic hydrogen by integrating through the disk and 
compare to observational limits. The line-of-sight hydrogen column density becomes: 



N(H) = 



323 



n 



13 



/ r \ 2A f r \ 
\ 117 AU/ + V 117 AU/ 



5.3 



1-1/2 



dr 



(2.1±5;g) x 10 20 



(10) 
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Integrating the distribution over the disk from 13 AU to 323 AU, we find that the total gas mass is: 

y323 

2it I I nn{r 1 z)r drdz 



H atoms 



13 
+ 17 
9 

+0.47 
0.08 



2 -2.9) x 10 bU atoms 



= (0.17™) Me. (11) 

The mass is within the limits from 2 1 cm observations (1.6 M m ; Freudling et al. 19951 ). but also within limits 
from dynamical arguments (0.4 iThebault fc Augereaul 120051 ). Thus we find that H I mass required to 
excit e the CO that we meas ure in absorption is close to the H I mass inferred in the solar composition case 
from IBrandeker et all ( 2004 ). 



DISCUSSION 



We observe CO absorption that is sub-thermally populated in the disk around j3 Pic. The observed 
properties of the absorption can be accounted for in two scenarios: one where H2O and CO are being 
photodesorbed off of icy grains, and an alternative scenario in which the CO that is coming off icy grains 
is excited by collisions with a residual gaseous disk of primordial hydrogen. The scenario involving H2O (or 
possibly another molecule) would result in a narrow region where the grains are being desorbed (<2.3 AU), 
while the scenario with residual, primordial hydrogen would result in a much larger region where the CO is 
being desorbed. 

The latter scenario requires that the CO and hydrogen are well mixed in the region where the CO is 
coming off the grains. The H I mass required to excite the CO is just at the 0.1 Mq required for ion- neutral 
collisions with hydrogen to explain the braking of the gas. Thus, it seems that while the excess carbon 
abundance is consistent with ion-ion collisions to be the breaking mechanism of the gas, the likely mass of 
the hydrogen in the disk is itself sufficient to brake the gas through ion-neutral collisions. The fact that the 
j3 Pic disk is so long-lived is not surprising given the fact that there is enough mass to brake the gas. 

If hydrogen is the dominant collision partner, we speculate that the metals and CO observed in the disk 
around /3 Pic are released from grains into a remaining reservoir of hydrog en gas. This may prov ide a possible 
explanation to the origin of the enhanced carbon abundance measured bv lRoberge et al] (|2006f ). Most of the 
metals (e.g. oxygen) may come from desorption of material throughout the disk as grains containing organic 
species with other metals desorb at the same rate as CO rtOberg et alJ l2009atfhh . However, if most of the 



oxygen is tied up in water while the carbon is in CH4 and CO, then the lower sublimation temperatures of 
those species compared to water (< 100 K versus 170 K for water) allow carbon-rich material to preferentially 
sublimate at the inner edge of the disk. Thus, carbon would be produced in excess relative to the other 
metals. 

The result presented here, especially the prediction for the m ass of hydrogen, can be tested with more 
sensitive measurements of hydrogen. Deeper observations than iFreudling et al.l (|1995) are necessary to 
test these predictions (an order of magnitude improvement would be sufficient to explore the estimated 
mass, while a factor of ^50 improvement could rule out the lower limit to 3cr) . The Square-Kilometer 
Aarray, estimated to be 50-100 times more sensitive than current radio telescopes (jCarilli fc Rawlingsll2004 : 
Dewdnev et al. 2009^ . should be able to test these predictions. 



Based in part on observations obtained at the Gemini Observatory, which is operated by the Association 
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of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of 
the Gemini partnership: the National Science Foundation (United States), the Science and Technology 
Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), Ministrio da Cincia e Tecnologia (Brazil) and SECYT (Argentina). 
The Phoenix infrared spectrograph was developed by the National Optical Astronomy Observatory. The 
Phoenix spectra were obtained as a part of programs GS-2002B-Q-43 and GS-2008A-C-7. Also based on 
observations obtained at the Infrared Telescope Facility, which is operated by the University of Hawaii under 
Cooperative Agreement no. NCC 5-538 with the National Aeronautics and Space Administration, Science 
Mission Directorate, Planetary Astronomy Program. M. R. T. acknowledges this work was performed under 
contract with the Jet Propulsion Laboratory (JPL) funded by NASA through the Michclson Fellowship 
Program. JPL is managed for NASA by the California Institute of Technology. 

Facilities: GeminkSouth (PHOENIX), IRTF (CSHELL). 
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Tabic 1. Log of j3 Pic Observations 



Date Instrument Spectral Range Integration Time S/N Lines covered 

(cm- 1 ) (m) 



2008 Mar 23 PHOENIX 2104 - 2113 16 90 v=l-0: P(8), P(9), v=2-l: P(l), P(2), P(3) 

2003 Jan 12 PHOENIX 2150 - 2160 8 100 v=l-0: R(l), R(2), R(3) 

2000 Aug 8 CSHELL 2146 - 2152 16 90 v=l-0: R(0), R(l) 



Table 2. Line Fluxes 



Line 




Vobs 


V ra d 


Equivalent Width 


N(CO) 


Flux 






(cm- 1 ) 


(cm" 1 ) 


(km s- 1 ) 


(xlO- 3 cm" 1 ) 


(xlO 14 cm- 2 ) 


(ergs cm- 2 s 




v=l-0 R(0) 


2147.08 


2146.95 


18 


4.9 ±0.8 


4.3 






v=l-0 R(l) 


2150.86 


2150.71 


21 


6.0 ±0.5 


9.0 






v=l-0 R(2) 


2154.59 


2154.43 


22 


4.4 ±0.5 


7.2 






v=l-0 R(3) 


2158.30 


2158.14 


22 


< 0.5 ±0.5 


< 2.3 






v=2-l P(l) 


2112.98 






< 0.5 




< 7.1 x 10- 


16 


v=2-l P(2) 


2109.14 






< 0.5 




< 7.1 x 10- 


16 


v=2-l P(3) 


2105.26 






< 0.5 




< 7.1 x 10- 


16 
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Fig. 1 — The spectrum of /3 Pic, including both the PHOENIX (< 4.651^m) and CHSHELL (> 4.651/zm) 
observations. The spectrum f3 Pic (solid line) and the telluric standard (dot-dashed line) are plotted in 
normalized units. The ratio of the spectra is also plotted and offset by 0.25 units. Areas where the trans- 
mittance is less than 50% are omitted. The locations of the CO transitions are shown for v=l-0. Absorption 
is observed from the low-J, v=l-0 levels. 
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Fig. 2. — The spectrum of (3 Pic. The spectrum f3 Pic (solid line) and the telluric standard, HR 5671 (dot- 
dashed line), have been normalized and overplotted. The ratio of the spectra is also plotted and offset by 0.25 
units. Areas where the transmittance is less than 50% are omitted. The locations of the CO transitions are 
shown for v=l-0 and v=2-l. No emission is observed, including the v=2-l P(l) line at ~ 4.733 /im. There 
appears to be a bump at the location of the v=l-0 P(9) feature, yet this is embedded in a combination of 
telluric features. Additionally, if the P(9) feature was a detection, the P(8) feature would also show emission, 
which is not observed. 
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13 AU 



323 AU 



Fig. 3. — Distribution of CO relative to other gaseous components in the /3 Pic syste m. The atomic 



component of the disk (H I and other metals) is observed to extend from 13 AU to 323 AU (jBrandeker et al 



2004 : Fernandez et al. 20061 ). The properties of the fundamental CO absorption reported here are consistent 
with an origin for this molecular component as CO photodesorbed from icy grains (gray objects) with an 
inner radius of 25 AU and an unknown outer radius. The CO is excited by collisions with the HI component 
of the disk. 
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Fig. 4. — CO excitation plot of the low-J ro-vibrational absorption lines. The data are shown as asterisks 
with error bars calculated from the ler variation of the continuum around the features, while the best-fitting 
model is overplotted as diamonds. The dashed line represents a temperature of 15 ± 2 K. The R(3) feature is 
an upper limit, and much lower than the LTE fit to the lower J lines, suggesting that J" = 3 is sub-thermally 
populated. The sub-thermal population allows the calculation of the density of the collision partner. 
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Fig. 5. — A one-dimensional chi-squared plot for the fitting of the equivalent widths using H2O as the 
collision partner. For each value of nH 2 o ; we calculated the rotational population of CO and the resultant 
equivalent widths of the v=0-l R(0), R(l), R(2), and R(3) lines. The reduced chi-square statistic for each 
run of our model is plotted, and the 68% confidence level corresponding to our la uncertainty is plotted as 
a solid line. We find that nH 2 o = (3.01^'q) x 10 3 cm~ 3 provides the best fit to the data. 
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Fig. 6. — A one-dimensional chi-squared plot for the fitting of the equivalent widths using hydrogen as the 
collision partner. For each value of njj, we calculated the rotational population of CO and the resultant 
equivalent widths of the v=0-l R(0), R(l), R(2), and R(3) lines. The reduced chi-square statistic for each 
run of our model is plotted, and the 68% confidence level corresponding to our la uncertainty is plotted as 
a solid line. We find that nn = (2.5j^' 2) X 10 5 cm -3 provides the best fit to the data. 
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Fig. 7.— Spectrum of CHSELL (> 4.651 //m) and PHOENIX (< 4.651 /im). The concatenation of the data 
from Phoenix and CSHELL reveals the R(0), R(l), and R(2) ro-vibrational CO absorption lines. The broad 
feature centered at 4.654 /im is the hydrogen line Pf (3. The spectrum has been modeled by calculating the 
excitation of the rotational levels in the ground vibrational state with hydrogen as the collision partner (that 
is, not assuming an LTE solution; see Section |4.3[) . The resultant absorption spectrum is plotted over the 
data (red dotted line). 



